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Abstract 

Background: Genetic epilepsy with febrile seizures plus (GEFS+) is a familial epilepsy syndrome with extremely 
variable expressivity. The aim of our study was to identify the responsible locus for GEFS+ syndrome in a 
consanguineous Tunisian family showing three affected members, by carrying out a genome-wide single 
nucleotide polymorphisms (SNPs) genotyping followed by a whole-exome sequencing. We hypothesized an 
autosomal recessive (AR) mode of inheritance. 

Results: Parametric linkage analysis and haplotype reconstruction identified a new unique identical by descent 
(IBD) interval of 527 kb, flanking by two microsatellite markers, 18GTchr22 and 15AGchr22b, on human 
chromosome 22ql3.31 with a maximum multipoint LOD score of 2.51. Our analysis was refined by the use of a set 
of microsatellite markers. We showed that one of them was homozygous for the same allele in all affected 
individuals and heterozygous in healthy members of this family. This microsatellite marker, we called 17AGchr22, 
is located in an intronic region of TBC1D22A gene, which encodes a GTPase activator activity. Whole-exome 
sequencing did not reveal any mutation on chromosome 22ql3.31 at the genome wide level. 

Conclusions: Our findings suggest that TBC1D22A is a new locus for GEFS+. 
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Background 

Epilepsy is one of the most common serious neurological 
disorders at worldwide level [1]. Environmental and gen- 
etic factors are known to play a role in its pathogenesis. 
Segregation studies suggest that most epileptic syndromes 
are complex disorders, but several monogenic forms have 
also been described. Genetic (formerly named generalized 
[2]) epilepsy with febrile seizures plus (GEFS+) is such a 
Mendelian inherited epileptic syndrome. This familial 
autosomal dominant (AD) epilepsy shows a wide range of 
phenotypes such as febrile seizures (PS), PS plus (PS+) - 
defined as PS persisting beyond the age of 6 - as well 
generalized and partial seizures [2,3]. In some families, 
most severe epileptic phenotypes have been described 
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such as the severe myoclonic epilepsy of infancy or 
Dravet syndrome [4-6] . 

Up to now, eight GEPS+ loci have been registered on 
OMIM database [7-16]. Three genes are known to be 
causative for GEPS+ phenotype: SCNIB (MIM#600235; 
GEPS+1) [7], SCNIA (MIM#182389; GEPS+2) [10] and 
GABRG2 (MIM#137164; GEPS+3) [11]. Pamilial forms of 
pure febrile seizures (PEB) have also been described. Segre- 
gation and twin studies of PS suggest a polygenic or multi- 
factorial mode of inheritance [17,18]. In most PEB families, 
phenotype follows an AD mode of inheritance [19]. To 
date, according to OMIM (Online Mendelian Inheritance 
in Man - http://www.ncbi.nlm.nih.gov/omim) database, 
eleven PEB loci have already been reported [20-30] and 
only four genes have been linked to PS related phenotypes: 
SCNIA (MIM#182389; PEB3A) [31], MASSl or GPR98 
(MIM# 602851; PEB4) [32], GABRG2 (MIM#137164; 
PEB8) [27] and CPA6 (MIM#609562; PEBll) [30]. 
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Finally, two other genes are considered as "suscepti- 
bility" genes to FEB syndrome and epilepsy: GABRD 
(MIM#137163; GEFS+5) [15] and SCN9A (MIM#603415; 
FEB3B/GEFS+7) [33]. 

In the present study, we reported the clinical and genetic 
studies of a GEFS+ consanguineous Tunisian family with 
an autosomal recessive (AR) mode of inheritance. We 
used homozygosity mapping, which is a method of choice 
for localizing genes responsible for AR diseases in consan- 
guineous families [34]. This approach allows identifying 
identical by descent (IBD) regions inherited from a 
common ancestor. IBD loci are homozygous in all the 
affected members but not in other relatives. In the present 
study we report a new locus on chromosome 22ql3.31 in 
which TBC1D22A gene (TBCl domain family, member 
22A, also known as C22orf4) is located. 

Results 

Clinical description of consanguineous Tunisian family 
with GEFS+ patients 

The familial pedigree is shown in Figure 1 and main 
clinical features of affected individuals are described 
in Table 1. The occurrence of generalized tonico 
clonic seizures (GTCS), absence seizure and FS+ led 



us to consider a GEFS+ syndrome in the present 
family. 

The proband (IV-4) is a 17-year-old boy, born at term 
after an uneventful pregnancy. At the age of 2, he 
suffered from a first complex FS. At the age of 6, a 
febrile GTCS and absence seizures were observed. Due 
to the continuation of FS, since the age of 5, he has been 
treated with valproic acid. Neurological examination 
revealed a mild mental retardation. Standard EEG tracings 
showed generalized spike waves. 

His first cousin, patient IV-3, is an 18-year-old boy, 
born at term after an uneventful pregnancy, from a con- 
sanguineous marriage between first degree cousins. In 
addition to FS+ which persisted until the age of 13, he 
experienced, at the age of 3, GTCS. Since the age of 8, 
absence seizures were also observed. Between the ages 
of 8 to 12, he was treated with valproic acid. Neuro- 
logical and EEG exams were normal. 

Patient IV-2, first cousin of the two previously de- 
scribed patients, is a 13-year-old girl, born at term after 
an uneventful pregnancy. She showed FS+ until the age 
of 10 and she suffered from GTCS since the age of 7. 
She has been treated by valproic acid since the age of 
5. Neurological examination and EEG tracings were 
normal. 
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Figure 1 Consanguineous Tunisian GEFS+family pedigree. Legend: A. Haplotype reconstruction for markers on chromosome 22q 13.31, 
which are ordered from centromere to telomere. Marl<ers from the lllumina array are in italic. The IBD haplotype shared by family members is 
shown in solid lines. Recombinaison events are shown by solid drawbar. Maternal and paternal chromosomes are designated by M and P, 
respectively. Each phenotype is describing by different colors. B. List of markers used for haplotype reconstruction with physical position. 
Markers from the lllumina array are in italic. 
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Table 1 Clinical characteristics of affected individuals 



Subject 


FS 




AFS 




EEG 


Neurological 


MR! 


Treatment/ 


(sex, age) 


Onset age/ 


n 


Onset age/ 


Type/n 




examination 




age 




remission, age 




remission, age 












IV-2 (F, 13y) 


4 y/yes, 1 0 y 


4 


7y/8y 


GTCS/3 


Normal 


Normal 


Normal 


VPA/7 


















y-until now 


IV-3 (M, 18y) 


2 y/yes, 1 3 y 


Numerous 


3y/7y 


GTCS/1 yearly 


Normal 


Normal 


Normal 


VPA/8 








8y/9y 


A/numerous 








y-12y 


IV-4 (M, 1 7 y) 


2 y/no 


Numerous 


6 y/1 1 y 


GTCS/5 


3 Hz generalized SW 


Mild mental 


Normal 


VPA/5 










A/numerous 




retardation 




y-until now 

















*AII individuals were agree to publisli clinical details. 

M: Male, F: Female, y: years, FS: Febrile Seizures, AFS: Afebrile Seizures, n: Seizures number, EEG: Electroencephalogram, GTC: Generalized Tonico Clonic Seizures, A: 
Absence, SW: Spike Wave, MRI: Magnetic Resonance Imaging, VPA: Valproic Acid. 



MRI scans were normal in all affected family members. 

Genetic linkage and haplotype analysis 

We assessed microsatellite markers (Additional file 1: 
Table SI) of four GEFS+ loci and seven FEB loci, that 
were known at that time. Haplotype reconstructions and 
analysis of recombination events allowed excluding all 
the GEFS+ and FEB loci analyzed in the present family. 

A genome wide genotyping was then performed and 
parametric linkage using the AR mode of inheritance 
identified a linked region of 10.6 Mb on chromosome 
22ql3.31 (Additional file 2: Figure FSl) flanked by the 
two SNPs rs3203726 and rs728592 with a maximum 
multipoint LOD score of 2.51 (Additional file 1: Table 
S5). On the same chromosome, positive LOD scores 
were also observed at six other loci (Additional file 1: 
Table S6). However, haplotype analysis showed that all 
these loci exhibited only one or two consecutive homo- 
zygous SNPs in every affected individual. Consequently, 
they are less likely to represent the disease locus. More- 
over, due to the number of obligate carriers of the 
disease (8) compared to affected ones (3), the AD model 
is not likely. 

The 10.6 Mb candidate region was further confirmed by 
using a set of microsatellite markers (Additional file 1: 
Table S2). For haplotype reconstructions, we used homo- 
zygosity mapping method, which allow localizing homozy- 
gous IBD genomic region(s) responsible for AR diseases in 
consanguineous families. We confirmed the existence of a 
common homozygous IBD haplotype shared by all the 
three patients. Therefore, the candidate region was re- 
duced to an interval of 527 kb, flanked by two microsatel- 
lite markers, 18GTchr22 and 15ACchr22b. This locus was 
homozygous in all three affected relatives and heterozy- 
gous in other family members (Figure lA). Within this 
homozygous IBD locus lie the 9 last exons and the 3'UTR 
of TBC1D22A which encodes the member 22A of Rab 
GTPase activator with a TBCl domain family (Figure 2). 
Two polymorphic markers, 17ACchr22 and rs2017931, 
are also displayed within the IBD interval (Figure lA). The 



most interesting one is the 17ACchr22 microsatellite 
marker due to its high variability. Its allele 4, correspond- 
ing to 16 AC repeats, was found homozygous in all af- 
fected relatives and heterozygous in all other unaffected 
family members. We estimated allele frequencies of 
17ACchr22 in a sample of 70 Tunisian controls. Ten 
alleles (corresponding to 13 to 22 AC repeats) were 
observed, and the most common alleles were allele 7 (19 
repeats) 25%, allele 8 (20 repeats) 22.1%, allele 5 (17 
repeats) 17.9%, and allele 6 (18 repeats) 10.7%. Finally, the 
IBD allele 4, had a frequency of 7.14% and was never 
found homozygous in the present Tunisian population 
(Additional file 1: Table S7). 

Array CGH analysis did not reveal any CNV at the 
genome-wide level in affected individual IV-3. 

Mutational analysis 

Mutation analysis was performed by whole-exome sequen- 
cing in patient IV-3 where 94.5% of the coding part of the 
RefSeq genes was covered at least 8x, which was sufficient 
to detect homozygous substitutions. This experiment 
revealed only three homozygous non-synonymous putative 
deleterious variations over the entire genome. Two of them 
were found on chromosome 22ql3.31, within the 10.6 Mb 
linked region: C.5159C > T (S1720L) in the TNRC6B (trinu- 
cleotide repeat containing 6B) gene and C.4433C > T 
(T1478M) (rs8141262) in the CACNAII (alpha II subunit 
of calcium voltage-dependant channel) gene. Genotyping 
these two missense variations, in all Tunisian family mem- 
bers, was done by HRM assay and confirmed by Sanger 
sequencing. These analyses showed that they are not linked 
to the disease phenotype, since the two other affected cous- 
ins (IV-2 and IV-4) were heterozygous for both variations. 
The third variation, found homozygous in patient IV-3, is 
located on chromosome lp36.11: C.604T > C (S202P) 
(rs6687605) in the LDLRAPl (low density lipoprotein 
receptor adaptor protein 1) gene. This polymorphism, 
genotyped by Sanger sequencing (Additional file 1: Table 
S3), was also found homozygous in his unaffected cousin 
IV- 1. Therefore, rs6687605 is not linked to GEFS+ 
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Figure 2 Structure of the TBC1D22A gene according to EnsembI (iittp://www.ensembl.org/index.iitml) withi localization of markers 
within and flanking the 527 Kb IBD linked locus. 
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phenotype in this family, according to the IBD assump- 
tion. Finally, we also performed Sanger sequencing for 
KCNJ4, since this gene codes for a ion channel protein. 
This family of genes are often found mutated in several 
epilepsies. Moreover, this gene is located within one of the 
6 chromosome 22 short homozygous regions (Additional 
file 1: Table S6), which are considered as not linked to 
the disease (see Genetic linkage and haplotype analysis 
section). No mutation was observed in this gene. 

Discussion 

The present study described clinical and molecular inves- 
tigations on GEFS+ syndrome in a Tunisian consanguin- 
eous family. The occurrence of both GTCS and absence in 
addition to FS+, and in absence of other types of seizures 
(e.g. myoclonic jerks), led us to consider that GEFS+ 
phenotype is likely within this family [2]. To our know- 
ledge, it would be the first description of a consanguineous 
GEFS+ Tunisian family with a putative AR mode of inher- 
itance, a transmission previously described in a Moroccan 
consanguineous family, showing FS and temporal lobe 
epilepsy phenotypes [30]. Moreover, all GEFS+ families, 
described in OMIM database, exclusively show an AD 
mode of inheritance [7-16]. Using homozygous linkage 
and whole-exome sequencing method, we mapped the 
pathogenic locus responsible for the disease within a short 
interval of chromosome 22ql3.31 in the present Tunisian 
consanguineous family. 

We first excluded by linkage analysis all known GEFS+ 
loci at that time, as well as those associated to FS pheno- 
types. Epileptic syndromes have also been observed in 
patients with structural variants [35,36], such as CNVs 
on chromosome 15ql3.3 are known to play an import- 
ant role in the genetic etiology of idiopathic generalized 
epilepsy [37]. We also excluded such genomic event in 
this family. 

The 22ql3.31 region, where we mapped the disease 
locus, has been implicated in a complex phenotype with 
epilepsy in 5-year-old boy with a 7.9 Mb de novo dele- 
tion of chromosome 22ql3.2-qter [38]. 

Based on the homozygosity of an IBD locus in patients 
only, refined mapping allowed reducing the disease locus 



to the 9 last exons and 3'UTR of TBC1D22A, and even to 
a shorter interval around the exon 12. However no muta- 
tion has been observed by whole-exome sequencing in 
any of these last exons of TBC1D22A. Homozygous alleles 
of the two polymorphic markers (17ACchr22 microsatel- 
lite and rs2017931), found in all patients, are likely to be 
non-functional and probably in linkage disequilibrium 
with the disease mutation in one of the TBC1D22A in- 
trons or in 3'UTR. The disease mutation could affect the 
splicing of TBC1D22A, Unfortunately, we were unable to 
test this hypothesis in absence of available RNA samples 
from the affected members of the family. 

The first crystal structure of TBC1D22A was done in 
2008 showing the diversity of human TBC domain fam- 
ily members [39]. Very recently, Shapshak et al. found 
that TBC1D22A gene had a differing gene expression 
profile across patients showing HIV associated dementia 
[40]. TBC1D22A gene is expressed in brain tissues with 
a percentage of 2.14% (http://smd.princeton.edu/cgi-bin/ 
source/sourceSearch). Moreover, GWAS suggested 
an association with human longevity and a SNP 
(rs5766691) located in TBC1D22A [41], within the 
527 kb IBD homozygous interval. 

Conclusions 

In conclusion, we reported a new locus on chromosome 
22ql3.31 in a consanguineous Tunisian family with a 
GEFS+ phenotype with an original AR mode of inherit- 
ance. We were not able to identify the disease mutation 
but restricted the linked interval to a part of TBC1D22A, 
which is expanded from exon 5 to the 3'UTR of the 
gene. Further studies will confirm whether this gene is 
also mutated in other families with GEFS+. 

This study was approved by the Research and Ethics 
Review Board of the Department of Neurology, University 
Hospital of Geneva and by the Ministry of Higher Educa- 
tion and Scientific and Technological Research of Tunisia. 

Methods 

Family ascertainment and phenotyping 

We recruited a consanguineous Tunisian family extended 
over four generations, comprising 3 affected and 19 non- 
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affected members, who all originated from Tunis (Tunisia) 
(Figure 1). The segregation of the trait suggests an AR 
mode of inheritance: all the affected relatives are on the 
last generation and one of them (IV-3) is born from first 
cousins (III-3 and III-4). The two other affected relatives 
were born from one parent with the same common ances- 
tors as III-3 and III-4. The two other parents (III-l and 
III-6) could descend from a putative common ancestor 
with their spouses because they all originated from a small 
village located in the North of Tunisia. 

All family members were agree to publish clinical details 
and they were clinically assessed at Charles Nicolle 
University Hospital (Tunis) by experienced epileptologists 
(Hela Khiari-Mrabet and Amel Mrabet). Information on 
FS and afebrile seizures, age at onset, duration, type and 
number of seizures, intellectual outcome, antiepileptic 
drug therapy, and seizures outcome were obtained from 
their parents and case notes when available. Electroen- 
cephalographic (EEC) and magnetic resonance imaging 
(MRI) scans were done for all affected family members. 
FS+ phenotype was defined as FS persisting beyond 
the age of 6 years [2]. Generalized tonic clonic seizures 
(GTCS) and absence seizures were defined according to 
the criteria established by the "International League 
against Epilepsy" [42-44]. 

Tunisian control group 

Seventy unrelated healthy Tunisian controls (mean age 
48.5 ± 16.5) were recruited at Charles Nicolle University 
Hospital (Tunis). These individuals did not have any per- 
sonal and familial history of seizures. 

Genomic DNA 

Oral and written informed consent was obtained from 
all participants or their legal representatives. Blood sam- 
ples were collected and genomic DNA was extracted 
from peripheral blood leucocytes by phenol/chloroform 
procedure [45]. 

Exclusion mapping 

Linkage analysis with a panel of 36 microsatellite 
markers spanning the FEB loci (FEBl, FEB2, FEB3, 
FEB4, FEB5, FEB6 and FEB11/ETL5) and GEFS+ loci 
(GEFS+1, GEFS+2, GEFS+3 and GEFS+4) known at that 
time were firstly explored (Table 2). Primer sequences 
are listed on Additional file 1: Table SI and PCR condi- 
tions are available as Additional file 1. We excluded all 
these previously linked loci by linkage analysis. 

Genetic linkage and haplotype analysis 

A genome-wide linkage study was performed by using the 
Illumina Gene Chip Linkage Infinium II HumanLinkage- 
24 Panel Breadchip® (Illumina, San Diego, CA), which 
allowed us to genotype 5913 SNPs. 



Table 2 Tested and excluded known FS and GEFS+ loci 



Locus 


Candidate 
gene 


Marker 


Genetic 
position 
(cM) 


Physical 
position 
(Mb) 


LOD 
Scores 


FEBl 




D8S553 


81.50 


67.10 


-9.31 






D8S1058 


90.10 


73.09 


-10.43 






D8S279 


90.20 


73.15 


-15.94 


FEB2 




D19S424 


10.97 


3.18 


-16.30 






D19S177 


20.75 


5.47 


-12.04 






D19S1034 


20.75 


6.06 


-8.99 






D19S406 


25.17 


7.33 


-13.27 






D19S76 


25.17 


7.56 


-9.23 


FEB4 


MASS] 


MASSllnt85 


99.30 


90.20 


-12.83 






D5S644 


104.76 


95.84 


-12.95 


FEB5 




D6S1620 


129.10 


129.99 


-5.31 






D6S472 


132.70 


132.58 


-9.52 


rbbo 


ll\ADA 1 


U 1 OO 1 1 JO 


o4./U 


1 U. 1 z 


- \ U.Uo 






IIVIrAzintj 


42.00 


1 2.01 


-1 6.96 






Dl 8S71 


42.80 


1 2.59 


-9.78 


r"iiiic 1 1 


jL/v/d 


u 1 yMz J 


jo./U 


40.1 9 


-11.14 






oLN 1 bint 1 


59.00 


40.21 


-1 1 .52 






P\ 1 nCODQ 

U 1 yboyo 


61 .40 


40.26 


-1 1 .43 


r'ccQ 1 ") 1 
uLro+z/ 

FEB3 


jLNI/V 
SCN2A 


UZjZddU 


1 75.50 


1 66.41 


-1 0.71 






D2S2345 


177.20 


168.43 


-6.59 






D2S2314 


188.90 


176.57 


-5.54 


GEFS+3 


GABRG2 


D5S1465 


162.00 


161.35 


-12.02 






GABRG2lntl 


162.50 


161.44 


-17.29 






GABRG2lnt5 


162.60 


161.48 


-17.29 






D5S2576 


1 62.63 


1 61 .51 


-1 6.1 7 






U 1 /OZ 1 D \ 


1 OZ./D 


1 ^; 1 QQ 
1 0 1 .oo 


- 1 j.yy 






UjMzz 


1 on 


1 Dz.uy 


-1 4.1 9 


r"IIIIC \ A 

uLrb+4 




p>oc 1 Q^n 
Uzb 1 ooU 


38.33 


1 7.36 


-3.68 






UZjdUd 


QQ Q"7 


\ y.zo 


-o.oy 






D2S2342 


40 47 


20 1 9 


-3 81 


FEBll/ 
ETL5 


CPA6 


D8S507 


75.00 


59.30 


-5.18 






D8S1812 


78.30 


60.85 


-5.47 






D8S1843 


78.80 


62.42 


-9.12 






D8S544 


81.00 


65.75 


-11.07 






D8S533 


81.50 


67.16 


-11.37 






D8S1775 


85.80 


68.99 


-11.99 



Two-point and multipoint logarithm of odds (LOD) 
scores were calculated with MERLIN 1.1 program, assum- 
ing an AR inheritance with complete penetrance, a disease 
allele frequency of 0.00001 and a phenocopy rate of 0. The 
haplotype reconstruction for family members was done 
manually, regardless of the individual affection status. 
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To confirm the region with a LOD score of 2.51 found 
by linkage analysis and in order to refine the mapping, we 
chose a set of 8 microsatellite markers located at the 
chromosome 22ql331 and distributed at average intervals 
of 0.2 Mb (Figure IB). According to the linked region 
defined by positive LOD score, we chose some microsatel- 
lite repeats using UCSC database (http://genome.ucsc. 
edu) by specif)dng "microsatellite" in "variations and 
repeats" part. Primer sequences (Additional file 1: Table 
S2) and PGR conditions are available on Additional file 1. 

Genome wide CNVs analysis 

We have also performed a copy number variations (CNVs) 
analysis by using the Agilent Human Genome compara- 
tive genomic hybridization (CGH) Microarray Kit 244 K 
(Agilent, Santa Clara, CA). The slide was scanned on an 
Agilent DNA microarray scanner. Data were obtained by 
Agilent Feature extraction software 9 and analyzed with 
Agilent CGH analytics 3.4 software, using the statistical 
algorithms zscores and ADM-2 with a sensitivity threshold 
of 2.5 and 6.0, respectively, and a moving average window 
of 0.2 Mb. Mapping data were analyzed on the human 
genome sequence using the NCBI database Build 35, 
Hgl7 (http://www.ncbi.nlm.nih.gov). 

Candidate gene mutational analysis 

Linkage results showed six other loci with positive LOD 
scores in the vicinity of the 527 kb IBD region. In one of 
them lies KCNJ4, which encodes a potassium inwardly- 
rectifying channel (subfamily J and member 4). According 
to its physical position and its functional role, this putative 
candidate gene for epilepsy, was sequenced in the present 
consanguineous family. We performed PGR for the two 
exons and splice junctions. We designed flanking exon 
primers from published sequences with the Primer3 
(version.0.4.0) online program (http://bioinfo.ut.ee/pri- 
mer3-0.4.0/primer3/input.htm). Primer sequences are listed 
on Additional file 1: Table S3 and PGR conditions are avail- 
able on Additional file 1. 

Exome sequencing 

Paired-end sequencing was performed using Illumina 
GAIIx/HiSeq 2000 instruments (Illumina, San Diego, CA) 
available at the Department of Genetic Medecine of the 
School of Medicine from Geneva. Sequence capture was 
performed using Agilent Sure Select Technology (Agilent, 
Santa Clara, CA) using either the All Exome kits or 
custom design corresponding to the linkage intervals and 
IBD regions. We typically obtain in excess of 30 million 
read pairs (60 million reads) of 76 nucleotides per GAIIx 
lane and 100 million reads pairs of 100 bp per HiSeq2000 
lane allowing multiplexing of 2-3 exomes per lane. 

Alignment of reads and call of single nucleotide variants 
and small indels were done using the latest version of 



MAQ/BWA [46] and PinDel [47] softwares. Only genetic 
variants with high-quality score were further investigated. 
Variant annotation was performed using the ANNOVAR 
package [48]. Generation of high-quality data necessary 
for reliable SNP detection across the targeted regions 
required 30-fold genomic coverage (equivalent to 15-fold 
coverage per haplotype). 

Each DNA variant was subsequently compared to Single 
Nucleotide Polymorphism database (dbSNP132, NCBI) and 
lOOOgenome project dataset. Only the unique and novel 
DNA variants were selected for additional filtering based 
on the nature of the mutation (gene structure location, nu- 
cleotide conservation, codon change, de novo occurrence) 
and familial segregation. Polyphen2 [49] and Sift [50] pro- 
grams were used to assess the pathogenecity of the filtered 
variants. If necessary, the final list of potential pathogenic 
variants were further refined by selecting genes func- 
tionally relevant to epilepsy. T1478M of CACNAII, 
S1720L of TNRC6B and S202P of LDLRAPl were further 
explored by classic Sanger sequencing and high resolution 
melt (HRM) assay using a Rotor-Gene 6000 instrument 
(Corbett Life Science, Australia). More details are given in 
Additional file 1 and in Additional file 1: Table S3 and 
Additional file 1: Table S4. 

Additional files 

f ^ 

Additional file 1: Table SI. Primer Sequer^ces of Microsatellite Markers 
for Exclusior^ Mapping. Table 52 Primers Sequences of Microsatellite 
Markers for Haplotype Analysis on 22q 13.31. Table S3 Primers Sequences 
for Sanger Sequencing. Table S4 Primers Sequences for HRM assay. 
Table S5 Lod Scores of linked locus on chromosome 22ql3.31. Table S6 
Loci on Chromosome 22ql3.31 with Positive LOD Score. Table S7 
17ACchr22 Microsatellite Genotype and Allele Distributions in the 
Tunisian Control Population. 

Additional file 2: Figure FS1. LOD scores were calculated using Merlin 
program. A new region was linked on the chromosome 22q with a 
maximum of LOD score of 2.51 . 
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